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The standard wave function approach for the treatment of neutrino oscillations fails in situations 
where quantum ensembles at a finite temperature with or without an interacting background plasma 
are encountered. As a first step to treat such phenomena in a novel way, we propose a unified 
approach to both adiabatic and non-adiabatic two-flavor oscillations in neutrino ensembles with 
finite temperature and generic (e.g. matter) potentials. Neglecting effects of ensemble decoherence 
for now we study the evolution of a neutrino ensemble governed by the associated Quantum Kinetic 
Equations, which apply to systems with finite temperature. The Quantum Kinetic Equations are 
solved formally using the Magnus expansion and it is shown that a convenient choice of the quantum 
mechanical picture (e.g. the interaction picture) reveals suitable parameters to characterize the 
physics of the underlying system (e.g. an effective oscillation length). It is understood that this 
method also provides a promising starting point for the treatment of the more general case in which 
decoherence is taken into account. 

PACS numbers: 13.15.+& 14.60.Pq, 14.60.St 



I. INTRODUCTION 



Ever since the conjecture that neutrinos have mass and thus might be subject to flavor oscillations [T| there has 
been a thriving interest in this very phenomenon which clearly indicates and gives rise to speculations as to how to 
C"| ■ describe physics beyond the electro-weak standard model. 

A convenient and established way to deal with neutrino flavor oscillations is to encode this effect in a Hamiltonian 
formulation in which the oscillatory behavior is captured in a Schrodingcr-like equation for a wave function in neutrino 
| flavor space. This formalism, in principle, applies to an arbitrary number of neutrino generations and is also capable 
, of incorporating medium effects on neutrino propagation such as coherent elastic forward scattering in, e.g., stellar 
matter [2|. It was soon realized that the Hamiltonian formalism for neutrino oscillations can be given a geometrical 
CTS , interpretation in N 2 — 1-dimensional flavor space for TV neutrino flavors Q. This approach to neutrino oscillations 
' sees equations of motion for a coherence vector [l5[ in that the Schrodinger-like equation of motion can be rephrased 
\ as a gyroscope equation, i.e. a formal equivalent to e.g. the precession of a magnetic moment in an external magnetic 
■ field. Besides its apparant usefulness when it comes to picture neutrino oscillations there is also a purely formal merit 
to the gyroscope-type equations in that they are introduced by means of decomposing the Hamiltonian in terms of 
the generators of the associated SU(N), e.g. the Pauli matrices for a two- flavor system. 

This decomposition procedure is also most convenient when the notion of a wave function is not suitable any 
more to describe the physics of neutrino oscillations. A typical situation in which the breakdown of the wave function 
formalism is expected are quantum ensembles with a finite temperature or neutrino ensembles with a finite temperature 
, and an interacting background plasma. The latter situation is encountered in the early universe prior to Big Bang 
Nuclcosythesis. The crucial point in such an environment is the breaking of coherence due to the small mean free path 
of neutrinos at high temperatures. The other important modification results from the fact that neutrino oscillations 
in the early universe can alter the lepton asymmetry, which in turn contributes to the refraction index of the primeval 
plasma rendering the equations of motion nonlinear. It is essentially for those two reasons that the wave function 
formalism must fail in describing neutrino oscillation phenomena now. The appropriate description is then given by 
the density matrix formalism. The density matrix of the neutrino ensemble obeys a von Neumann equation and the 
different contributions to the effective Hamiltonian are given by collisions (non- forward scattering), with particles 
from the background medium, which introduce decoherence. On the other hand coherent oscillations are governed by 
a matter-dependent effective Hamiltonian which comes about via coherent forward scattering processes of neutrinos 
off the background particles [1, H[ ■ 
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The dynamics of the neutrino ensemble are determined by the Quantum Kinetic Equations (QKE) which present 
a generalization of the Pauli-Boltzmann equations. The former evolve quantum amplitudes as is indispensable if 
a consistent description of particle oscillation phenomena, which are inherently nonclassical, is sought. The Pauli- 
Boltzmann approach on the other hand evolves probabilities rather than amplitudes; this procedure is essentially 
classical since quantum mechanics only enters the problem when it comes to calculating cross sections for the various 
possible reaction channels. The resultant Quantum Rate Equations (QRE) are inappropriate when neutrino oscilla- 
tions occur. Thus, in order to obtain the QKE the full density matrix for all particles in the plasma is evolved forward 
in time by means of the S matrix and tracing over all degrees of freedom other than the neutrinos under consideration 
yields the equation of motion for the system's density matrix, the QKE, which do reduce to QRE in the appropriate 
limit [H]. The variable of interest in the QKE is the one-body reduced momentum-dependent density operator, which 
is conveniently decomposed in terms of the generators of the associated SU(N). 

In the forthcomig analysis we account for a brief motivation on how to obtain the QKE from the density matrix 
formalism in the case of coherent forward-scattering (which also dominates the bulk of the studies to follow) and how 
to relate the solutions to the QKE, i.e. the coherence vector, to physical observables in section [Til Moreover, we 
introduce the Magnus expansion, which allows for an analytic, yet approximative, solution to the underlying QKE. A 
short discussion of the convergence properties of the Magnus expansion proves useful to single out suitable quantities 
to describe the physics at hand. In section IIIII we unfold a systematic way to develop an adiabatic perturbation 
theory starting from the QKE for a two-flavor neutrino ensemble with generic potentials. We show that the Magnus 
expansion not only allows to analytically solve the QKE in a perturbative way, but can also serve to motivate the 
definition of physical quantities, such as an effective mixing angle or an adiabaticity parameter. In order to isolate 
a convenient perturbation parameter we perform different changes of the quantum mechanical picture (e.g. into the 
interaction picture) for the QKE. Hence, a suitable succession of bases changes can improve the convergence properties 
of the expansion. We solve the QKE for the coherence vector to first order in the Magnus expansion. The perturbation 
parameter in the adiabatic perturbation theory can then be used to identify a suitable expansion parameter for the 
non-adiabatic case. Once this parameter is identified the non-adiabatic perturbation theory can be treated on the 
same ground as the adiabatic one, i.e. singling out an appropriate expansion parameter by changing the quantum 
mechanical picture and solving the QKE in the resulting representation by means of the Magnus expansion. This 
non-adiabatic perturbation theory is developed in section HVl In section [V] we comment on the integrals appearing in 
both the adiabatic and non-adiabatic perturbation theories and show how the latter lead to the correct limits such as 
vacuum neutrino oscillations and the slab model approximation to non-adiabatic transitions in the MSW framework. 
We sketch how to extend the developed perturbation theory to higher orders in the Magnus expansion. Moreover, we 
elaborate on possible extensions of the theory introduced to scenarios including dccohcrence in the ensemble due to 
collision processes with particles from a background plasma. It turns out that the continuation of the perturbation 
theory is somewhat nontrivial and deserves a more careful study. 



The purpose of our analysis is to describe the evolution of a two-flavor neutrino ensemble with generic potentials 
at a finite temperature. To this end we commence our considerations by a close inspection of the underlying QKE, 
which in this context can be readily derived from the density matrix formalism for neutrino oscillations. We remark 
that the latter has to be used to correctly deal with quantum systems at finite temperature as well as in situations 
where loss of coherence becomes important. The wave function formalism ceases to provide an appropriate handle 
for such systems. However, it is not until section IY CI that wc discuss an extension of the formalism to be discussed 
shortly to systems with a background plasma, hence decohering collisions, present. 

The density matrix p(p 1 1) obeys a von Neumann equation 



where p is the neutrino four-momentum, H a generic Hamiltonian for the system. Since we are dealing with a two- 
flavor system we can now decompose both the density matrix and the Hamiltonian in terms of the generators of the 
associated SU(2), namely the Pauli matrices er\ This yields 



II. NOTATION AND MATHEMATICAL TOOL BOX 



f>(p, t) 



i[H(p,t),p(p,t)] 



(1) 




(2) 



(3) 



where repeated indices are to be summed over and the obvious identifications Pq = Trp, Vb = TriJ, P = Tr (per), 
V = Tr (Ha) and a = (a x , <r v ,<r z ) have been made. Here P is the so-called coherence vector Q. Making use of this 
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notation it is straightforward to recast the von Neumann equations according to 

Pi = [-V e ilk V k ] P l . (4) 



Hence we can identify 




Su = -V e Uk V k or S = V \ V z -V x (5) 



as the evolution matrix of the neutrino ensemble [l6| . Put another way, the QKE can now be written [TtJ as 

±P(t) = S(t)P(t). (6) 
The entries of the effective potential vector V — (V x , V y: V z ) can also be readily obtained via 

A rrr 

V V X = p = 2ReH 12 = -—sin 2^o, (7) 

2p 

V V y = 2Im# 12 = 0, (8) 

Am 2 

V V Z = A = (ff u -H 22 ) = —cos 28 + V a , (9) 

2p 

where Hij are the elements of the Hamiltonian H , Am 2 is the mass-squared difference of the two neutrino states, 9q 
is the associated vacuum mixing angle between the two flavors, which we shall denote as v a and z/f, for defmiteness. 
V a is the difference of potential terms affecting v a and respectively. The last equality has been obtained using the 
2x2 neutrino oscillation Hamiltonian in flavor space. 

In the coherence vector description the expectation values of the generators of the associated SU{2) are promoted 
to observables of interest. All information about the system can thus, in principle, be extracted from a solution to 
Eq.© for P. Some comments related to this issue are in order: In the one-particle interpretation the diagonal entries 
of the density matrix simply give the probability to find the system in one or the other state, i.e. 

prob (i/„ i/„) = -P [l + P z ], (10) 
prob (!/„-»> i^) = ip [l-P.]- (11) 

In the ensemble interpretation of the density matrix the diagonal entries give the relative number densities for the 
different neutrino flavors normalized to the equilibrium Fcrmi-Dirac number distribution at zero chemical potential fi 
according to 

N a (p) = ±P [l + P z )N E( Z(p,0), (12) 
N b (p) = ±P [l-P z )N E( Z(p,0), (13) 

NEQ (p^) = i^r^' (14) 

Z7T x — (- e T 

where T is the temperature of the ensemble. Those relations are also easily inverted to yield a physical meaning of 
both 

N a + N b 

Po = N m ■ ( 15 ) 

N a N b 

2 N a + N b - (W) 

Hence Pq is connected to conservation of probability and, in a broader context, also lepton number. It is important 
to note that oscillations merely swap neutrinos from one flavor to another so that Pq does not evolve in time, unless 
repopulation effects from some background plasma have to be taken into account as is the case, e.g., in the early 
universe. On the other hand P z parametrizes the asymmetry of the system that is the excess of v a over v b . The latter 
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fact also motivates the way of speaking in which P x , P y are called coherences encoding the amount of decoherence in 
the system. Therefore the evolution of P z is of special interest in most applications. 

Our next concern is to provide mathematical means to approximately solve the differential equations (J5]) ■ In order to 
do so we rewrite the QKE supplied with an initial condition as 

^P(t)=S(t)P(t), P*=P(to) (17) 

and notice that we are dealing with a non-autonomous set of linear differential equations. If the matrix S did not 
depend on time, the corresponding differential equation could be readily solved by taking the matrix exponential 
exp[S(t — to)] and writing P(t) = exp[S(t — to)]P l . It is then fair to ask whether the solution to Eq.(fT7|) can always 
be written as an exponential via 

P(t) =exp[fi(t,t )] P l - (18) 

A method for finding such a true exponential solution [l8| has indeed been established under the label Magnus 
expansion 0, H| . The Magnus operator Q = In S satisfies its own differential equation which in turn is solved by 
a series expansion fl = X^j^Li ^n(t) , where each Magnus approximant £l n {t) is small in an appropriate sense. The 
smallness of each Magnus approximant then determines the convergence properties of the expansion. A general 
theorem states that for a differential equation (fT7j) defined in a Hilbert space with a bounded operator S the series 
converges in the interval t 6 [to , t c [ such that 

dr ||5(r)|| < 7T, (19) 

where 1 1 . 1 1 is a matrix norm. As can be inferred from this condition the convergence properties of the Magnus expan- 
sion can be improved by means of a change of the quantum mechanical picture and thus improving the convergence 
properties of the expansion can be rephrased as minimizing the matrix norm of the evolution matrix. 
The first two Magnus approximants assume a form 



f drS(r), (20) 

Jt 

fia(t) = \ I dh I' dt 2 [S(h),S(t 2 )] (21) 

1 Jto Jta 

and various methods have been worked out to calculate higher order terms @- All such terms contain nested 
commutators of S evaluated at different times from which it is readily seen that the Magnus expansion gives the exact 
result (the matrix exponential for a time- independent matrix S) already in first order, if the matrix S commutes with 
the matrix obtained by integrating S over a certain time interval. Note that each approximant Q n adopts the same 
properties as the matrix S. If 5* was (anti-)Hermitian, so would be the exponentials of all approximants; this means 
that unitarity of the solution will be preserved in each order of the perturbation expansion seperately. 

The convenience of the Magnus expansion is now apparent: In lieu of a calculation of eigenvalues and eigenvectors 
of a matrix it poses the problem of calculating a matrix exponential. The latter problem is usually solvable for 3 x 3 
matrices. Note, however, that the dimension of the evolution matrix corresponds to the number of generators of the 
associated SU(N), i.e. for an A-dimensional system the dimension of the evolution matrix is N 2 — 1. So for the 
case of three interacting neutrinos one already has to choose between diagonalizing a 8 x 8 matrix or calculating its 
exponential. Both problems are intricate. 

Another nice feature of the Magnus expansion is that it provides a clear description of how to improve the approxi- 
mation by going to higher order terms. It also reproduces exact solutions to the QKE in the appropriate limits as will 
be seen later, but is also apparent at this stage of our analysis by inspecting the form of the Magnus approximants. 
Moreover, the condition for its convergence appears as an integral condition on the norm of the evolution matrix. 
This condition can be used as a cross check for whether changing from one quantum mechanical picture to another 
improves convergence properties of the expansion or not. We will invoke this criterion in the upcoming analysis and 
as it turns out, it can even supply some physical insight. 
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III. ADIABATIC PERTURBATION THEORY 




Having established QKE of the form 

|p(t) = S(t)P(t), S(t) = I A(t) -/}(*) I (22) 

it is easy to calculate the matrix norm (l9j according to 

||S||! = Tr(StS)=2u4. (23) 

It is straightforward to show that the effective oscillation length of the system is indeed given by 

2tt 

T e 



— = UeS =y/#+p. (24) 



At this point of our analysis it might seem academic to calculate the matrix norm of the evolution equation under 
consideration, but it will be seen shortly that a comparison between matrix norms in different quantum mechanical 
pictures can provide physical insight into the nature of the neutrino ensemble at hand. Note, moreover, that we 
treat both j3 and A as time-dependent quantities. This generic assumption provides greater freedom when it comes 
to adapting the QKE to early universe applications. To this end we notice that j3 scales as p _1 in momentum and, 
in an expanding universe, this momentum is rcdshifted and thus depends on time. We shall sketch how to treat such 
situations later on. 

Furthermore, in order to get a grasp on how this oscillation length can be understood physically we transform the 
QKE to a basis which resembles the commonly encountered mass eigenbasis in the MSW framework. To do so we 
notice that there are only non- vanishing contributions to the effective potential vector's x and z component, namely 
V x and V z . Thus it is only sensible to consider a generic time-dependent rotation in the ccz-plane by an angle 9(t) as 

cos9(i) sin9(t)\ 

P(t) = R(t)Q(t) with R[@(t)} = ( 1 (25) 

•sin9(t) cos 6(f)/ 

where Q is the coherence vector in the new co-rotating frame and R(t) is the time-dependent rotation matrix. The 
QKE in the new basis appear as 

B ( -A cos 6-/3 sin 6 

— Q(t) = S Q (t)Q(t), Sq= Acos6 + /?sinG -/3cos9 + Asin9 I . (26) 



dt 



^ P cos 8 - A sin 9 



Since we introduced 9 as a generic time-dependent mixing angle we are endued with its explicit definition according 
to our needs. It is readily seen that the (Sq) 23 and (Sq) 32 elements of the evolution matrix can be eliminated by an 
appropriate choice of the mixing angle 9(i). The advantage of this choice is the geometrical interpretation: in the 
Q-picture the motion of the coherence vector is confined to the xy-plane, if there was not the additional perturbation 
by the time derivative of the effective angle, which introduces a non-zero z-componcnt to the problem and forces the 
motion to exit the xj/-plane as the ensemble evolves. The smaller the change of the effective mixing with time, the 
smaller the urge of the coherence vector to exit the xy-plane. Therefore we fix the effective mixing angle to be 

cos 6ft) = — A(t) sin9(t) = -= ^ (27) 

The effective mixing angle reveals that mixing becomes maximal (9 = tt/2) if the condition 

A(* res ) = (28) 

is satisfied for the so introduced resonant time t res . A vanishing A(£), i.e. maximal effective mixing, hence coincides 
with the existence of a resonance in neutrino conversions, which can also be equivalently rephrased for a resonant 
temperature T rcs , depending on the application one has in mind. 
We now recast the evolution matrix in the Q-picture as 



Sq = Weff . (29) 
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Consider the matrix norm of this evolution matrix in the new quantum mechanical picture 



\ S q\\1 = 2w cff 



1 



i de N 

W c fF dt 



(30) 



At first glance the above transformation seems to worsen the convergence properties due to the appearance of the 
additional 

7= -77 (31) 

w c ff dt 

term. However, if this very term is sufficiently small, 7 <C 1, the convergence will only be marginally altered. Moreover, 
the smallness condition can be understood physically as well: the characteristic time scale of the system under study 
is T S y S = 1/ Weil] whereas the characteristic time scale of the interaction can be identified as Tint = (dG/dt)- 1 . Hence 
the parameter 7 simply compares the characteristic time scale of the system to the characteristic time scale of the 
interaction, stating that a small 7 can be paraphrased as the system's time scale being much smaller than the 
interaction's time scale. Put another way, the interaction is adiabatic. The parameter 7 is thus readily interpreted 
and henceforth referred to as the adiabaticity parameter for the system. 
The evolution matrix of the system thus reads 




S Q = w eff . (32) 



Before we move on with our analysis it is just to briefly comment on the definition of the adiabaticity parameter and 
the effective mixing angle. The effective mixing angle defined above is essentially the expression encountered when it 
comes to the usual MSW framework of matter-affected neutrino oscillations. Note, however, that the latter typically 
features a sin 20 instead of sin as defined here. In order to streamline notation we will nonetheless still omit this 
factor of two. Moreover, we alert the reader that defining I/7 as the adibaticity parameter is also quite common in 
the literature. However, the physics is not altered by this convention. Also when comparing our analysis to other 
work it is important to notice that in the instead of 20 convention the adiabaticity parameter lacks a factor of two 
as well. We will analyze the adiabaticity parameter further in section HVl 

It is by now established that 7 can serve as a small perturbation parameter in the adiabatic regime of neutrino 
conversions. It thus feels harmonious to struggle through just another transformation which is introduced to isolate 
the perturbation parameter 7 in a convenient way and such that fast convergence of the expansion to come is assured. 
The prescription is as follows 



$(t) = U(t)X(t), where ^U(t) = S%(t)U(t), U(t ) = 1, (33) 



having decomposed the evolution matrix according to 

Sq = S% + S^ q = (oj cS I + ( I (34) 

in self-obvious notation. The subsidiary evolution equation for U(t) is also readily solved to give 

(cos (D ff — sinw c ff 0\ nt 
sinweff cosweff and 0J c ft(t) = / dr w c ff(r). (35) 

I J Jt o 

The QKE in disguise are recognized to be 
d 



— JLU e ff COS UJ e ff ' 

7w ff sii 

,7w off cosw off -7w Gff sina) off 



— X{t) = S x (t)X(t), Sx{t)=\ 7w cff sin£ cff | (36) 



and calculating the matrix norm yields 

||5x|||= 2 7 2 Wc 2 ff - (37) 



7 



It is evident now that the small parameter in the adiabatic regime, namely 7, has been isolated and hence good 
convergence of the sought-after perturbation theory can be expected. This comforts us to seek a perturbative expansion 
in this basis (which is but an interaction picture for the Q-basis). 

The considerations unfolded in this section have seen two linear transformations R(t), U (t) from the original P-basis 
to the Q- and X-basis. The reason for those transformations is twofold: On the one hand changing the basis for the 
QKE discloses the physics of the system we are dealing with and on the other hand it seems advisable to find a basis 
for the QKE in which an approximate solution gives accurate results. For convenience we shall now recapitulate the 
meaning of the transformations introduced so far. 

The first transformation (P — ^ Q) is inherently physical. It gives a recipe how to establish the concept of a mass 
eigenbasis in the coherence vector description of neutrino oscillations. The effective mixing angle defined in this way 
differs from the effective mixing angle encountered in the common MSW formalism by a conventional factor of two. 
In this quantum mechanical picture a clear path of approaching the resonance in neutrino oscillations is unveiled. 
A resonant conversion of neutrino flavors is encountered for X(t) = 0. Moreover, the transformation to the matter 
eigenbasis sees the introduction of an effective mixing angle, which is a harbinger for the adiabaticity parameter 7 
subsequently layed open. Adiabatic neutrino conversion occurs for 7 <C 1, when the time scale of the system is much 
smaller than the time scale of the interaction. The mathematical benefit of this transformation is that we get a grasp 
on the convergence properties of the approximation we want to employ and we can furnish it with physical meaning. 
The convergence properties of the expansion get worse as the amount of adiabaticity violation increases. A fact which 
later will be useful to construct a non-adiabatic perturbation theory. Also the change to the matter eigenbasis suggests 
that the adiabaticity parameter 7 should be the appropriate small quantity to expand in. 

The second transformation (Q — > X) is convenient from a mathematical point of view. It removes an exactly 
intcgrable part of the evolution matrix and thus the matrix norm for the latter is directly proportional to the small 
expansion parameter 7. This truly renders 7 into the sought-after perturbation parameter and we take comfort that 
the envisaged expansion converges fast. Put another way, already the first approximant should provide good guidance 
for the exact solution. 

Note eventually that no attempt for solving the QKE has been made so far. We have merely changed the quantum 
mechanical pictures to unfold the underlying physics. The paradigm of our analysis is that a careful treatment of the 
QKE, i.e. a succession of different changes of the quantum mechanical pictures already allows to extract important 
information about the system under consideration without ever explicitly solving the QKE. 

In the A-basis the QKE are solved to first order in the Magnus expansion by 



exp 



dr S x (t) 



to 



X(t ). (38) 



The formal solution for the coherence vector P{t) to first order in the Magnus expansion is thus obtained as 



P^\t)=R{t)U{t) 



exp 



/ dr Sx(t) 

J t 



R-\t )P(t Q ). (39) 



In order to streamline notation we write the terms contained in the matrix exponential as 

J s (t) = / dr 7w cff sinw ff, (40) 

J t 

J c (t) = / dr 7w off cos lj oS (41) 

Jto 

as well as 

I J| ee JJTTJI (42) 
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and the resultant expression for the coherence vector assumes a form 

cos6(£) sin0(i)\ /cosw c ff — sinews 0^ 

P (1) (i) = ( 1 J sinweff cos(I> off 

sin6(t) cos6 / I 1 



1 



J 2 + J c 2 cos|J| —J C J S (—1 + cos | J|) -J c | J| 2 sine | J| N 

|2 



X— 2 -J c J s (-l + cos|J|) J 2 + J 2 cos|J| J s \J\smc\J\ 
J c |J| 2 sinc|,/| -J s |J| 2 sinc|J| |J| 2 cos|J| 

^cos 6 — sin 8 ^ 

1 | P\ (43) 
^sinOo cosOo 

where additionally 6 (to) = ©o and sinca; = was defined. This expression, as cumbersome as it may look at a 
first glance, presents an analytic, yet perturbative, solution to the QKE as given in Eq. (|22| for a generic potential, 
i.e. for a generic time (or equivalently temperature in early universe applications) dependence of both j3 and A, as 
long as the transition can be considered adiabatic (7 <C 1). Note also that oscillation probabilities in the one particle 
interpretation can lie extracted from this formal solution. Moreover, oscillating contributions to this very probability 
can be studied since there is no inherent averaging over rapidly oscillating contributions as is usually considered in 
the derivation of the oscillation probability in the MSW framework. Still, to fully appreciate this result a thorough 
discussion of various limiting cases, such as the adiabatic limit, is called for [To| . We postpone this endeavor until 
section [V] 



IV. NON- ADIABATIC PERTURBATION THEORY 



Before we proceed by developing a non-adiabatic perturbation theory on the same grounds as the foregoing adiabatic 
perturbation theory it is instructive to briefly reconsider the adiabaticity parameter as defined in Eq. pip and rephrase 
it in a way that allows for an understanding of the notion of adiabaticity in terms of the parameters (3 and A. It is 
straightforward to show that 

. . /3A d /3 

However, physically the adiabaticity parameter at the neutrino conversion resonance (A = 0) is of foremost interest. 
We find 



1 dA 



(45) 



The adiabaticity parameter depends on the shape of the matter profile dA/di as is expected from the MSW framework. 
Large variations of the matter profile at resonance are clearly disfavored for the sought-after perturbation expansion 
to work. Besides this contribution also the term 1//3 2 is familiar. It simply states that the larger (3, the better the 
expansion works. Put another way, if the notion of adiabaticity as put forward in our analysis is adopted, a small (3 
at resonance is incompatible with an adiabatic perturbation expansion to some extent. A regime with small /? thus 
indicates non-adiabatic transitions; we will also refer to this regime as the sudden regime henceforth. 

Recalling the QKE according to Eq.© it is obvious that /3 itself can be adopted as a small perturbation parameter 
if a non-adiabatic perturbation theory is desired. We split the evolution matrix S as 

-A(t) 0\ /0 
S(t) = S x (t) + S/}(t) = \ X(t) + -j9(i) I (46) 

0/ \0 P{t) 

in obvious notation. The S\ subsystem of this evolution equation can be integrated exactly and hence we impose the 
following change of the quantum mechanical picture 

P(t) = V(t)Y(t) with ^V(t) = S x (t)V(t). (47) 
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The subsidiary evolution equation again is solved to give 

''cos A — sin A 0^ 
V(t) = | sin A cos A 

1; 



and 



A(i) 



dr A(t), 



(48) 



to 



rephrasing the QKE as 



-Y(t) = Sy(t)Y(t), 



S Y (t) 






y/3 sin A /? cos A 




(49) 



Note also, that the new quantum mechanical picture is just the interaction picture. The matrix norm reveals isolation 
of the small perturbation parameter: 



>y | If 



= 20 1 



(50) 



Solving the QKE for Y{t) will see the time-integrated evolution matrix Sy since the first order Magnus expansion 
gives the coherence vector as 



P (1) (t) = V(t) 



cxp 



dr Sy{t) 



to 



P(to) 



(51) 



and it is thus sensible to define the following integrals to streamline notation 



K s (t) 
K c (t) 



dr {3 sin A, 



dr /3cosA 



as well as 



\K\ = y/K*+K* 



(52) 
(53) 

(54) 



in order to mimic the notation introduced above for the adiabatic perturbation theory. The coherence vector to first 
order in the Magnus expansion is calculated to be 



'cos A —sin A 0^ 
sin A cos A | x 

1; 



1 



K 2 + K 2 cos\K\ K c K s (-l + cos | K \) -K s \K\ 2 sine \K\ 



\K 



a \K C K S (-1 + cos \K\) K 2 + K 2 cos\K\ -K c \K\ l sine \K\ \P l 



K, \K\ 2 sine \K\ 



(55) 



K r K sine K 



K cos \K 



This is the analytic perturbative solution to the QKE in Eq. (f22|) with generic potential and time dependence for f3 and 
A as long as the evolution can be considered non-adiabatic, which is equivalent to saying that /3 is a small quantity 
one can expand in. Again, this result can only be fully appreciated once the associated limiting cases are recovered. 
We will see to this in the next section Hill. 



V. PERTURBATION THEORY INGREDIENTS AND LIMITING CASES 



We understand that our approach is a generic solution to the QKE in Eq. f)22f) without making any explicit reference 
to the physics it can be applied to. Hence the integrals J c / S (t), K c / S (t) explicitly depend on the time dependence of 
both A, i.e. the potential term V a , as well as (5 and have to be evaluated in each application separately. However, 
certain general statements can be made already due to the fact that neutrino conversions reveal a resonance at A = 0. 

In any case, however, it is still necessary to demonstrate that the Magnus expansion does give exact results in the 
various physical limits. 
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A. The integrals: J c / S (t) and K c / a (t) 

The structure of the integrals J c / S (t) and K c / S (t), for the adiabatic and non-adiabatic case, respectively, at a first 
glance, is similar: the integrand is the expansion parameter (juj c s for the adiabatic case; /3 for the non-adiabatic case) 
multiplied by an oscillatory function. On second thought, however, there is a crucial difference; the oscillatory term 
in K c / S has a stationary phase (dX/dt = at resonance), whereas J c / S does not. 

Let us therefore evaluate K c / S by means of the stationary phase method. Two main assumptions are needed to 
apply the stationary phase method: the oscillatory behavior of the integrand must be rapid enough to suppress all 
large contributions to the integral which might come from /3 so that the latter can simply be evaluated at resonance. 
The other requirement is that the resonance happens in a small region around i res ; put another way, the smallness of 
the aforementioned region is determined by whether the substitution to — > — oo and t — > oo is justified in this region 
or not. Supposing that these two requirements are met, we find 



hence 



\K S \ 

\Kc\ 

\K\ 



0. 



2tt 

^res 
^rcs 



1/2 



1/2 



(56) 
(57) 

(58) 



Note, that in the non-adiabatic perturbation theory the reciprocal value of the adiabaticity parameter at resonance 
is a small quantity. 

Consider the integrals J c / S now. The first step that comes to mind here is integration by parts. We get 



Js(t) 
Jc(t) 



-7 cos cjcff 



7smw c ff 



dr 



to 



d7 
d7 



COS W e ff, 



, d7 . _ 
dr — sin cj e g . 
dr 



(59) 
(60) 



If the variation of 7 in the interval [to,t] is sufficiently mild the main contribution to the integrals is expected to come 
from the first term on the right hand side. 

Given these arguments the integrals K c / S , J c i s reveal a common trademark. Both integrals turn out to be small in 
the following sense: K c / S is proportional to the inverse of 7 ros which, in a non-adiabatic perturbation theory is a large 
quantity; likewise in the adiabatic perturbation theory 7 is the small quantity to expand in and again the integrals 
J c / S turn out to be proportional to 7. 



B. Limiting cases 

We begin our consideration by studying appropriate limits for the adiabatic perturbation theory first. 

1. Adiabatic perturbation theory: The vacuum limit 

This limit to Eq. (|43[) is probably the most intuitive one, if we confine ourselves to the one particle interpretation. 
We understand that for an exactly solvable system the first order Magnus term should already give the exact 
result, hence implying P^\t) = P(t). Let us examine how this works out here: Firstly, we discard the potential 
term V a . This gives A — > — Ap- cos2$o and j3 — ^p- sin2#o- Hence, for consistency we must take the adiabatic 
limit of 7 — > 0, which immediately implies J c / S — > 0. Moreover, a) e g — > ut if we define the common oscillation 
frequency oj = 2 " and set to = (as no resonance time exists, the choice of to is arbitrary). Finally, it follows 
directly from the definition of the effective mixing angles that 6(t) — > 29o for all times. All this reduces the 
coherence vector to 

(cos 2 2# cos ujt + sin 2 26 cos 2# sin u)t sin 4# sin 2 ^ \ 

cos 26» sin wt coswt - sin 26> sin ut \ P\ (61) 

sin 40 o sin 2 ^ - sin 28 Q sin cut cos 2 26 + sin 2 26 cos wt J 
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Suppose we start with a v a flavor such that P l = (0, 0, 1). The probability to find the neutrino in the same/the 
other state after time t is then using Eqs. (|10HTT]) as well as Eq. (|6"Tj) given by 

prob(^ a ->■ v a ) = 1 - sin 2 26*o sin 2 (62) 

prob(^ a — > v ) = sin 2 26 a sin 2 ^- (63) 

which is just the common oscillation probability. Note, moreover that this result was obtained solely using the 
truncated Magnus expansion as given above and that it accounts for probability conservation. Put another way, 
unitarity is guaranteed by means of the expansion itself and does not have to be imposed by hand. 

We next turn our attention to the non-adiabatic case. There exist (at least) two interesting limits. 

1. Non-adiabatic perturbation theory: The sudden limit 

The sudden limit of taking /? — )► renders the QKE (|2"2"|) into formally exactly solvable differential equations 
such that the Magnus expansion should give an exact result. The limit /3 — > enforces w c g = |A| and hence 
cos O = 1, sinG = 0, which in turn implies 7 — > 0. The coherence vector assumes a form 

''cos A — sin A 0\ 

P{t) = I sin A cos A O P*. (64) 
1/ 

The coherences of the ensemble are oscillating as a function of time (the ensemble is incoherent) and the flavor 
is frozen to its initial value. Put another way, in physical situations in which the evolution of the ensemble 
happens in a way that with increasing time also /3 increases, the unfreezing of the ensemble can be studied using 
non-adiabatic perturbation theory since it treats j3 as a small perturbation. We will point out in section IV CI 
that this is typically the case in early universe applications. There is, however, a twist when it comes to early 
universe applications in that such systems are typically collision-dominated at high temperatures and thus the 
notion of adiabaticity is expected to be modified due to the presence of collisions. Put another way, a small (3 
in early universe environments augmented by the presence of decohering collisions might as well allow for an 
adiabatic perturbation theory (see section IV CI for some more details). 

2. Non-adiabatic perturbation theory: The slab model limit 
Suppose that we start the evolution of the neutrino ensemble from a purely v a state P % = (0, 0, 1). We obtain 
for the coherence vector 

f (—K s cos A + K c sin A) sine \K\ 
P (1) (*) = ( -(K c cos A + K s sin A) sine \K\ \ ■ (65) 
cos\K\ 

The flavor oscillation probability is then written as 

prob(v a -)• v a ) = cos 2 i \K\ . (66) 

Suppose furthermore that the situation as described above to estimate the K-type integrals holds, i.e. the reso- 
nance in neutrino conversions happens in a narrow time interval centered around i rGS . Applying the stationary 
phase approximation we then obtain 



.2 1 



prob(f a ->• v a ) = cos 4/ (67) 

V 2 7r CS 

This result is the oscillation probability for the slab model as outlined in Ref. [l2j]. So the relevant limit is 
respected in this situation. 



C. Higher order corrections and applications 



As has been seen in the previous sections the Magnus expansion can be easily extended to higher orders by summing 
the associated approximants according to f2(t) = £7i (t) + Q,2(t) + In order to get a grasp on how this prescription 
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unfolds we calculate the 2nd order Magnus approximant to be 

/ J 

n 2 (t) =\-J o o | , 

where J is given by 




Jad = r / dti / dt 2 7(ii)7(t 2 )wcff(ii)wcff(i2)sin[a) off (t2) - w c ff(^i)] (69) 



for the adiabatic case and 



t /.ti 

to 



Jnad = £ / 



for the non-adiabatic case respectively. The calculations are performed in the A-picture for adiabatic transitions 
and in the K-picture for non-adiabatic transitions. Two things are readily inferred: the second order approximant is 
indeed 0(j 2 ) and O(0 1 ) for adiabatic and non-adiabatic corrections respectively as is expected. Moreover it is seen 
that the 2nd order populates the (23) and (32) entries of the Magnus operator fl. 

We shall now shortly elaborate on the complications which arise when ensemble decoherence is to be taken into 
account. This typically happens in early universe applications in which the time evolution of neutrinos is governed by 
three distinct physical processes: Firstly, the expansion of the universe. Secondly, coherent oscillations governed by 
a matter-dependent effective Hamiltonian which results from coherent forward scattering processes of neutrinos off 
the background particles. Thirdly, scattering processes with the background plasma of elementary particles. These 
collisions, or non-forward scattering processes, with particles from the background medium typically introduce deco- 
herence effects into the neutrino ensemble. In our analysis in this paper we have neglected the ensemble decoherence 
due to non-forward scattering. 

The epoch of foremost interest in studying neutrino oscillations in the early universe is the one between muon 
decoupling at T ~ m M ~ 100 McV and neutrino decoupling, i.e. prior to Big Bang Nucleosynthesis (BBN), at about 
T ~ 1 MeV, since during this time the initial conditions for Nucleosynthesis, the electron neutrino abundance, are 
set (2^1 which then directly influence the neutron-to- proton ratio at the onset of BBN via /3-processes p + e~ ^ n + v e . 
The primordial plasma during this epoch thus consists of electrons, positrons, neutrinos and antineutrinos. 

The density matrix p for the system of interacting and oscillating neutrinos encodes ratios of number density 
distributions and hence the expansion of the universe does not directly contribute to the time evolution of the density 
matrix. However, the momenta of the particles are redshiftcd and the equilibrium number distributions N E Q(p, 0) 
thus depend on time through this rcdshifting. 

So far most attention has been paid to two neutrino systems and specifically to active-sterile flavor oscillations 
[l3| . The latter are interesting since active-sterile oscillations would populate the additional sterile species and thus 
contribute significant additional energy density, which in turn would trigger an accelerated expansion of the universe 
and hence lead to a higher weak freeze-out temperature. This again would alter the neutron-to-proton ratio at the 
onset of BBN. The coherent part of matter-affected active-sterile oscillations splits in two contributions. One is just 
the leading order density-dependent contribution (the Mikheev-Smirnov-Wolfenstein 0] (MSW) part). This part 
is only temperature-dependent indirectly via the cosmological rcdshifting of fermion number density. The second 



contribution comes from leading order finite temperature gauge boson effects [F|, which cannot be neglected at the 
temperatures considered here. 

The loss of coherence is due to neutrino collisions with the background medium. The decoherence (or damping) 
function for this process in thermal equilibrium turns out to be proportional to the total collision rate for the neutrino 
with momentum p under consideration. 

The epoch of interest can now be decomposed into three distinct domains: At high temperatures finite temperature 
gauge boson effects dominate. Rcpopulation effects from the background plasma can be neglected since at high 
temperatures thermal equilibrium for all relevant species is rapidly established. At intermediate temperatures lepton 
number production starts and the forward scattering contribution comes into play as a small perturbation. Finally 
prior to the onset of BBN, at low temperatures, collisional effects and finite temperature gauge boson contributions 
cease to be important and coherent neutrino oscillations are the dominant process. 

In each of the aforementioned temperature domains the QKE should be solved to determine the evolution of the 
neutrino ensemble. It is obvious that decoherence and rcpopulation effects modify the underlying QKE and complicate 
their analysis by introducing new physical scales in the system. The early universe framework thus deserves a more 
careful treatment which is beyond the scope of this work. 
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VI. CONCLUSIONS 



We have developed a new formalism to solve the Quantum Kinetic Equations governing the evolution of a two- flavor 
neutrino ensemble with finite temperature subject to generic potentials. To this end we use the Magnus expansion 
which has virtues when it comes to pcrturbatively solve coupled non- autonomous differential equations such as the 
QKE considered here: the expansion is unitary to each order of the approximation and its convergence properties are 
determined by the matrix norm of the underlying evolution matrix, which depends on the physical parameters of the 
system. The perturbation ansatz proposed maps the problem of diagonalizing a Hamiltonian with generic potentials 
(and hence possible time dependences) onto calculating the exponential of a matrix and the associated integrals as 
its emergent ingredients. 

Contrary to previous approaches the formalism presented here describes both adiabatic and non-adiabatic flavor 
conversions in a quantum neutrino ensemble at finite temperatur and arbitrary (e.g. matter) potentials. In the 
adiabatic case a suitable succession of bases changes is employed to define an effective oscillation frequency, mixing 
angle and adiabaticity parameter for the ensemble. The latter is then used as a perturbation parameter. In the 
non-adiabatic case the results obtained for the adiabaticity parameter are used to identify its inverse at the resonance 
as a suitable expansion parameter. We give explicit expressions for the coherence vector to first order perturbation 
theory for both adiabatic and non-adiabatic conversions in the ensemble. We show that our approach is compatible 
with the appropriate limiting cases for both adiabatic and non-adiabatic conversions and elaborate on how to tackle 
the integrals emergent from the perturbation theory. 

We understand that our analysis presents a promising base for further investigations concerning the inclusion of 
decoherence in the ensemble and the resulting applications in early universe scenarios, e.g. the generation of lepton 
number in active-sterile oscillations prior to BBN. Moreover, the approach may be useful for the analysis of long- 
standing oscillation anomalies such as MiniBooNE or even MSW-like scenarios in astrophysical environments. 
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